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ABSTRACT 


Laminar/turbulent  transition  is  studied  in  a  straight  channel  with  1  Hz 
imposed  bulk  flow  unsteadiness  for  Reynolds  numbers  of  1 103  to  2715  and 
Strouhal  numbers  of  0.06045  to  0.02455.  Channel  aspect  ratio  is  40  to  1  with 
1.27  cm  height,  50.8  cm  width  and  427  m  length.  Observations,  videos,  and 
photographs  of  smoke  patterns  show  different  subcritical  transition  events 
including  Tollmien-Schlicting  waves,  vortex-array  type  motion  evidenced  by 
smoke  swirls  and  ribbon  like  patterns,  turbulent  spots  and  fully  turbulent  flow. 
Fluctuating  intensity  magnitudes  relative  to  velocities  from  imposed 
unsteadiness  are  determined  from  phase-averaged  velocity  traces.  Of 
particular  interest,  is  whether  the  imposed  unsteadiness  advances  or  delays 
transition  events.  With  the  conditions  that  were  investigated,  transition 
events  begin  to  occur  at  lower  Reynolds  numbers  and  extend  over  a  wider 
range  of  Reynolds  numbers  compared  to  flows  with  no  imposed  unsteadiness. 
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L  INTRODUCTION 


A.  BACKGROUND 

Many  typical  flows  found  in  nature  and  in  aerodynamic  applications  and 
almost  all  flows  associated  with  turbomachinery  involve  freestream 
unsteadiness.  Most  of  these  flows,  especially  those  of  technological  interest 
undergo  transition  from  a  laminar  state  to  a  turbulent  condition.  The  accurate 
prediction  of  the  location  and  extent  of  transition  and  the  accompanying 
changes  of  important  flow  properties  is  vital  to  the  design  of  aircraft  and  most 
energy  transformation  machines  involving  fluid  flow. 

Laminar-turbulent  transition  results  in  significant  increases  in  local  wall 
shear  stress  and  convective  heat  transfer  rates.  This  becomes  evident  by  first 
considering  laminar  flows,  where  transport  of  heat  and  momentum  normal  to 
streamlines  occurs  only  as  a  result  of  molecular  diffusioa  This  is  usually  an 
advantage  since  tangential  forces  at  solid  surfaces  (drag)  and  heat  flux  normal 
to  surfaces  are  relatively  low.  Occasionally,  laminar  flows  are  undesirable  since 
they  tend  to  separate  from  surfaces  more  easily  than  turbulent  flows.  In  most 
applications,  laminar  flows  are  unstable  and  transitional  or  fully  turbulent  flows 
exist  In  the  latter,  higher  rates  of  diffusion  across  mean  streamlines  due  to 
turbulent  transport  cause  significantly  higher  drag  and  heat  flux  at  surfaces. 

The  parameter  generally  used  to  describe  the  onset  of  transition  is  the 
Reynolds  number.  In  channel  flows,  linear  stability  theory  predicts  that  the 
laminar  flow  becomes  unstable  to  small  amplitude  disturbances  at  a  critical 
Reynolds  number  of  5772  Here  the  Reynolds  number  is  based  on  mean 
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centerline  velocity  and  channel  half-width.  With  large  amplitude  disturbances, 
however,  the  transitional  Reynolds  number  for  channel  flow  may  be  as  low  as 
1000- 1500.  The  value  of  Reynolds  number  at  the  onset  of  transition  depends 
strongly  on  the  disturbances  existing  in  the  flow  and  conditions  which  prevail 
in  the  entry  region  of  the  channel. 

Also,  the  sequence  of  events  which  occur  during  transition  are  strongly 
dependent  on  upstream  history,  especially  laminar  boundary  layer  character 
and  development  There  is  a  growing  collection  of  evidence  which  also 
suggests  that  turbulent  shear  layer  development  is  more  strongly  dependent 
on  its  upstream  history  than  previously  believed.  This  would  mean:  (1)  that 
transition  affects  turbulent  boundary  layer  structure  downstream,  and  (2)  that 
turbulent  boundary  layer  results  obtained  under  equilibrium  conditions  (Le. 
streamwise  self-similarity)  may  not  adequately  represent  flow  behavior  near 
surfaces  of  operating  turbines  and  airfoils.  Non-equilibrium  shear  layer 
development  near  turbine  passage  surfaces  and  near  airfoils  is  due  to  varying 
surface  curvature,  varying  pressure  gradients,  and  freestream  unsteadiness. 
Unsteadiness  is  particularly  important  in  this  regard  since  it  affects  shear 
properties  at  a  given  location  in  two  ways:  ( 1)  by  its  local  direct  influence  and 
(2)  from  its  affect  on  shear  layers  upstream,  whose  history  affects  downstream 
development  The  computational  work  of  Singer,  Ferziger  and  Reed  [Ref  1J 
indicates  that  moderately  low  frequency  oscillations  in  the  channel  can  initiate 
nonlinear  effects  which  trigger  transition  on  flow  regimes  where  the  steady 
flow  remains  laminar.  These  happen  even  though  linear  theory  indicates  that 
the  flow  is  stable.  Experimental  clarification  of  these  points  is  needed,  along 
with  the  verification  of  the  sequence  of  events  occurring  during  the  transition 


process  at  different  Strouhal  numbers.  The  view  of  the  overall  transition 
process  is  especially  valuable  to  the  gas  turbine  research  community. 

Unsteadiness  is  additionally  important  since  boundary  layers  which  exist 
near  turbulent  blades  are  significantly  different  from  equilibrium  boundary 
studied  in  many  laboratories.  A  number  of  recent  studies  indicate  the 
inadequacy  of  results  obtained  from  equilibrium  boundary  layers.  Gillis  and 
Johnston's  (Ret  2]  measurements  in  boundary  layers  downstream  of  a  convex 
curved  surface  and  Ligrani's  [Ref  31  results  in  boundary  layers  downstream  of 
artificial  thickening  devices  showed  that  many  structural  changes  persist  long 
distances  downstream  of  perturbations.  This  indicates  that  boundary  layer 
behavior  may  depend  on  upstream  history  to  a  greater  extent  than  previously 
believed  This  is  important  regarding  unsteadiness  since  it  alters  boundary 
layers  locally  as  well  as  the  history  of  boundary  layers  upstream. 

Few  experiments  have  been  performed  to  study  the  stability  or  transition 
of  oscillating  flow  in  a  plane  channel.  Most  known  information  is  from  a 
handful  of  analytic  studies  and  the  full  Navier-Stokes  simulation  undertaken  at 
Stanford 

Regarding  the  latter,  the  work  by  Singer,  Reed  and  Ferziger  [Ref.  1] 
provides  fascinating  information  as  to  the  eifects  of  unsteadiness  on  transition 
For  plane  channel  flow  linear  theory  has  been  verified  at  a  Reynolds  number  of 
5000.  With  mean  flow  frequency  in  the  range  of  the  non-dimensional 
frequency  of  the  least  stable  TS  wave  (Tollmien-Schlichting),  they  find  that 
the  mechanisms  driving  the  transition  are  similar  to  those  described  by 
Herbert's  secondary  instability  theory  [Ref.  4  and  5)  and  observed 
experimentally  by  Nishioka  et  aL  [Ref.  6  and  7],  and  Kozlow  and  Ramazanov 


[Ret  8J.  In  the  plane  channel  flow  these  mechanisms  have  been  studied 
extensively  in  previous  computational  investigations.  The  addition  of  mean 
flow  modulation  causes  significant  variations  in  the  amplitude  of  the  TS  wave 
during  the  period,  hence  the  dominant  flow  structures  depend  on  the 
Strouhal  number  and  the  magnitude  of  the  modulation. 

In  a  straight  channel  without  imposed  unsteadiness  the  process  of 
transition  begins  with  TS  waves  experiencing  mild  growth  at  sufficiently  high 
Reynolds  numbers.  When  these  waves  grow  to  an  amplitude  of  about  0.3 
percent  of  the  free  stream  velocity,  three-dimensional  waves  become 
unstable.  Transition  often  occurs  at  Reynolds  numbers  below  that  at  which  TS 
waves  become  unstable.  These  waves  cause  Lambda  shaped  vortices  to  lift 
from  the  near  wall  region  and  ■burst*  ejecting  low  speed  fluid  from  close  to 
the  wall  into  the  faster  moving  fluid  away  from  the  wall.  These  events  then 
lead  to  a  rapid  breakdown  of  the  laminar  boundary  layer  into  a  turbulent  one. 
The  process  can  be  modified  by  a  number  cf  different  factors.  Freestream 
turbulence  can  lead  to  non-sinusoidal  oscillations  in  the  nominally  laminar 
boundary  layer.  Local  disturbances,  such  as  roughness  elements,  may  lead  to 
structures  known  as  turbulent  spots.  These  patches  of  turbulent-!  ike  flow 
form  a  characteristic  arrowhead  shape  and  spread  in  the  downstream 
direction.  Distributed  roughness,  acoustic  disturbances,  surface  waviness  and 
other  factors  all  serve  to  complicate  the  basic  process.  There  also  exist 
equilibrium  states  where  two  dimensional  waves  (which  are  similar  in  form  to 
TS  waves)  with  a  given  amplitude  will  neither  grow  nor  decay.  The  initial 
disturbances  in  the  inlet  of  the  channel  can  evolve  into  these  or  similar 
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structures,  after  which  the  strong  secondary  instability  of  three-dimensional 
waves  can  take-  place. 

B.  OBJECTIVES 

This  research  is  a  follow-on  to  the  work  undertaken  by  Longest  (Ret  9), 
and  is  believed  to  be  the  first  experimental  (quantitative)  study,  known  to  the 
author,  on  the  effects  of  imposed  mean  flow  unsteadiness  on  transition  in 
straight  channel  flow.  The  objective  is  to  investigate  the  effects  of  imposed, 
periodic,  and  deterministic  unsteadiness  on  the  transition  process  from 
laminar  to  turbulent  flow.  Answers  to  the  following  questions  are  sought 

1.)  What  parameters  govern  the  unsteady  flow? 

2)  What  events  occur  in  the  flow  during  the  onset  and  development  of 
transition? 

3. )  What  is  the  interaction  between  the  organized  induced  unsteadiness 
and  transition  (e.g.,  development  of  two-dimensional  Tollmien-Schlichting 
waves  into  three-dimensional  instabilities)? 

4. )  At  what  range  of  frequencies  can  one  expect  interactions  between 
the  organized  imposed  unsteadiness  and  transition  phenomena? 

Flow  measurements  (using  hot-wire  anemometry)  and  flow  visualization 
are  used  to  analyze  the  effect  of  oscillating  the  bulk  flow  in  otherwise  steady 
plane  Poiseuille  flow.  An  attempt  is  made  to  determine  whether  the  type  of 
unsteadiness  studied  (sinusoidal  at  one  and  two  hertz)  contributes  to  the  early 
onset  of  transition,  d“lay  of  transition,  or  has  no  effect  Sinusoidal  oscillation 
was  chosen  to  induce  a  deterministic  periodic  variation  of  velocity.  The  results 
are  compared  with  the  results  of  Singer's  [Ref.  1)  numerical  simulation,  which 


which  suggests  that  such  oscillations  provide  a  stabilizing  effect  at  all  but  very 
low  frequencies. 

C.  ORGANIZATION 

Section  U  follows  this  introduction  and  describes  the  experimental 
facilities.  In  section  III  experimental  procedures  are  discussed,  including 
calibrations,  channel  validation,  flow  measurements  and  flow  visualizatioa 
The  phase  averaging  technique  is  also  detailed  here.  A  discussion  of  the 
results  follows  in  section  IV.  Conclusions  are  presented  in  section  V. 
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EL  EXPERJM  1 MTAL  FACILl  i  lfiS 


A.  CHANNEL 

The  40  to  1  aspect  ratio  straight  channel,  shown  schematically  in  Figure 
1,  is  located  in  the  laboratories  of  the  Department  of  Mechanical  Engineering. 
The  details  of  the  design  and  construction  are  discussed  by  Ligrani  and 
Subramanian  (Ref  lOj.  Briefly,  the  channel  is  made  of  6.35  mm  thick 
plexiglass  with  a  straight  test  section  4.27  m  in  length,  with  inside 
dimensions  of  1.27  cm  in  height  and  50.8  cm  in  width.  The  top  and  bottom 
walls  of  the  channel  are  supported  by  ribs  and  cross  beams  along  its  length. 
The  side  wall6  are  removable  in  order  to  gain  access  to  the  inside  of  the 
channel.  The  inlet  section  consists  of  a  honeycomb  assembly,  3  frames  with 
screens  and  a  20: 1  contraction  ratio  nozzle.  Two  layers  of  cheesecloth  are 
placed  at  the  inlet  lip  for  filtering  the  air  entering  the  channel.  At  the  exit  of 
the  4.27  m  long  test  section  are  three  10.16  cm  long  frames  with  screens  for 
housing  the  unsteady  device  ( the  device  used  to  create  unsteadiness  in  the 
flow)  and  a  second  honeycomb.  A  two  dimensional  diffuser  45.72  cm  long 
with  a  3  degree  total  angle  is  located  just  upstream  of  the  plenum  chamber. 
The  exit  plenum  chamber  has  inside  dimensions  of  60.96  cm  by  60.96  cm.  A 
3.81  cm  orifice  plate  assembly  between  the  plenum  chamber  and  the  blower  is 
provided  to  meter  the  flow  through  the  channel.  The  exit  of  the  plenum 
chamber  is  connected  to  the  suction  side  of  a  5  HP.  blower  with  a  long  50.8 
cm  diameter  pipe.  A  valve  is  provided  between  the  plenum  and  the  blower  to 
regulate  the  flow.  Bulk  flow  velocities  up  to  15  m/s  are  achieved  in  the 
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channel.  In  our  experiments,  bulk  flow  velocities  in  the  range  1.0-45  m/s 
are  used. 

The  channel  is  designed  with  a  sufficient  length  such  that  transition 
occurs  after  the  laminar  velocity  profile  (parabolic)  is  fully  developed  (self 
similar).  The  initial  flow  development  length  is  estimated  (Ref  1  l:p.  186]  to 
be  within  about  25  per  cent  of  the  test  section  length  for  the  range  of  speeds 
used  in  the  present  tests. 

B.  UNSTEADY  DEVICE 

The  most  important  aspect  in  the  design  of  any  unsteady  device  for  a 
transition  study  is  to  ensure  that  the  device  produces  a  deterministic  and 
periodic  unsteadiness  without  adding  other  disturbances.  To  achieve  this  in 
an  open  circuit  induction  tunnel,  the  unsteady  device  is  best  located  at  the 
exit  of  the  channel.  This  way,  problems  caused  from  the  wakes  and  the  flow 
disturbances  resulting  when  the  unsteady  device  is  upstream  of  the  test 
section  are  avoided. 

The  design  is  based  on  the  design  used  by  Miller  and  Fejer  [Ref.  12]  and  is 
shown  in  Figure  2  Since  the  depth  of  the  channel  is  only  127  mm.  a  single 
rotating  vane  driven  through  a  spur  gear  train  by  a  DC  stepper  motor  is 
adequate  to  introduce  the  required  unsteadiness  in  the  flow.  The  rotating 
vane  causes  varying  resistances  to  the  flow  in  the  test  section.  The  frequency 
of  the  imposed  unsteadiness  is  controlled  by  changing  the  vane  rpm  and  the 
amplitude  is  controlled  by  using  different  vane  widths.  It  is  also  important  to 
match  the  flow  speed  with  the  vane  speed.  A  vane  8.7  mm  wide  was  used  in 
this  study. 
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The  vane  in  Figure  2  is  made  of  a  3.2  mm  thick  brass  strip  with  rounded 
edges,  and  spans  the  entire  width  of  the  channel  The  vane  is  supported  at 
the  ends  by  a  3.2  mm  diameter  shaft  and  bushings  that  are  fitted  to  the  side 
walls  of  the  frame.  Three  additional  intermediate  spanwise  struts  are 
provided  to  increase  the  rigidity  of  the  vane.  One  end  of  the  shaft  is  extended 
to  accommodate  a  48  TPI  (Threads  Per  Inch)  spur  gear.  The  driver  spur  gear 
has  12  TPI  and  is  mounted  on  the  Superior  Electric,  M092-FD310  Stepper 
Motor.  To  achieve  a  fine  step  control  a  Modulynx  MITAS  PMS085-D050 
Drive  with  Modulynx  MITAS  PMS085-C2AR  Drive  Controller  is  used.  This 
combination  facilitates  motor  shaft  speeds  up  to  100  revolutions  per  second, 
with  a  constant  200  steps  per  motor  shaft  revolution.  Note  that,  with  a  1:4 
gearing  ratio  and  two  cycles  of  imposed  flow  oscillations  per  vane  rotation,  one 
motor  revolution  corresponds  to  one  half  cycle  of  the  imposed  flow  oscillation 
by  the  vane.  The  Drive  Controller  can  be  programmed  for  the  required  motor 
speed  (rps).  The  Modulynx  PMS085-D050  Drive  was  also  used  to  provide  a 
trigger  signal  to  initiate  data  acquisition.  This  trigger  also  enabled  the 
correlation  of  data  taken  with  the  position  of  the  vane  at  the  time  that  the  data 
was  recorded 

C  FLOW  MEASUREMENTS 

The  equipment  used  for  flow  measurement  includes  a  digital  manometer, 
hot-wire  probe,  hot-wire  bridge,  signal  conditioner,  digital  multimeter,  high 
speed  data  acquisition  system,  oscilloscope  and  a  microcomputer.  These  are 
described  in  detail  below. 
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L  Bulk  Velocity 

The  bulk  flow  velocity  through  the  channel  is  determined  from  the 
pressure  drop  across  the  3.81  cm  orifice  plate  which  is  installed  in  the  piping 
between  the  blower  and  the  exit  plenum  chamber.  This  pressure  drop  is 
measured  and  displayed  by  a  Validyne  Model  PS  309  digital  manometer 
(inches  of  water). 

2.  Instantaneous  Velocity 

A  constant  temperature  hot-wire  anemometer  is  used  to  measure 
instantaneous  velocities  of  the  flow.  Briefly,  the  principles  of  hot-wire 
anemometry  (constant  temperature  mode)  are  as  follows. 

The  flow  sensing  element  (hot-wire  probe)  consists  of  a  short  length  of 
wire,  connected  to  two  supports,  which  forms  the  leg  of  a  Wheatstone  bridge. 
This  bridge  maintains  the  wire  at  elevated  temperature  and  is  initially 
balanced  prior  to  exposing  the  probe  to  the  flow.  This  balancing  resistance, 
referred  to  as  "cold"  resistance  is  multiplied  by  an  "overheat  ratio".  Once  the 
bridge  is  set  to  this  value,  and  as  the  probe  is  exposed  to  the  flow,  the  bridge 
is  placed  in  operation  and  a  small  current  (60  mA)  passes  through  the  circuit 
Wire  temperature  increases  because  of  resistive  heating  and  the  wire 
resistance  increases  until  the  bridge  is  rebalanced.  The  governing  relationship 
is: 

R  =  Ro(l  +a <T  -  To)) 

where  R  is  the  resistance  at  temperature  T,  Ro  is  the  resistance  at  the 
reference  temperature  To.  and  a  is  the  temperature  coefficient  of  resistance. 


Flow  past  the  probe  causes  the  wire  to  cool  (by  convection)  and  the 
bridge  becomes  unbalanced  as  the  resistance  decreases.  The  bridge  is 
rebalanced  by  a  feedback  loop  and  amplifier.  Response  is  rapid  enough  to 
maintain  wire  temperature  and  resistance  nearly  constant  As  such,  the 
variations  in  the  bridge  voltage  are  directly  proportional  to  the  current 
variations,  and  flow  velocity  can  be  calculated  as  a  function  of  the  bridge 
voltage. 

Now  the  different  components  used  to  measure  instantaneous  velocity 
are  described. 

a.  Hot-wire  Probe  and  Probe  Positioning 

A  Dantec  Electronics  Inc.  P51  probe,  mounted  such  that  the  wire  is 
horizontal  and  normal  to  the  streamwise  flow  direction,  is  used.  The  wire 
sensor  is  5.0  microns  in  diameter  and  approximately  1.25  mm  long.  The 
probe  is  mounted  through  the  side  wall  of  the  channel  and  is  held  fixed  by  a 
rotatable  lever  arm  which  allows  placing  of  the  probe  at  various  y/d  positions 
ranging  from  0.2  to  03.  A  scale  graduated  in  terms  of  y/d  is  mounted  adjacent 
to  the  lever  arm  for  probe  positioning  to  an  accuracy  of  approximately  0.5  mm. 
Longitudinal  and  transverse  position  of  the  probe  for  this  study  were  46.12  cm 
and  15.64  cm.  measured  from  the  near  flange  at  the  end  of  the  test  section 
and  from  the  inside  of  the  left  channel  wall  (looking  downstream), 
respectively.  By  rearranging  the  sidewalls  the  probe  can  be  positioned  at 
other  streamwise  locations  for  future  studies. 

b  Hot-wire  Bridge 

The  hot-wire  is  operated  with  a  DISA  55  M10  constant  temperature 
bridge.  A  five  meter  coaxial  cable  is  used  to  connect  the  hot-wire  to  the 


bridge.  For  the  purposes  of  this  study  the  cold  resistance  was  multiplied  by  an 
overheat  ratio  of  1.8.  The  DC  voltage  from  the  bridge  is  measured  with  a 
Hewlett-Packard  3466A  digital  multimeter. 
a  Signal  Conditioner 

The  hot-wire  signal  is  filtered  and  amplified  by  a  DISA  56W20  signal 
conditioner.  The  high-pass  filter  is  set  to  0.1  Hz  in  all  cases  to  remove  the  DC 
signal.  This  however,  did  not  alter  the  effects  of  1  and  2  Hz  imposed  flow 
unsteadiness.  The  low-pass  filter  is  set  at  10  KHz  to  filter  out  noise  existing 
above  this  frequency.  The  hot-wire  signal  is  amplified  using  a  gain  of  10  to 
achieve  an  acceptable  signal  to  noise  ratio.  Output  from  the  signal  conditioner 
is  sent  to  the  high  speed  data  acquisition  system  and  to  a  B  &  K  Precision  10 
MHz  Oscilloscope  used  to  signal  monitor  the  signaL 
d  Highspeed  Data  Acquisition  System 

A  Hewlett-Packard  6944A  Multiprogrammer  is  used  in  conjunction 
with  a  Hewlett-Packard  series  9000,  Model  310  microcomputer  to  acquire 
data  A  buffered  69759A  A/D  card  configuration  is  used  with  the 
Multiprogrammer  to  allow  continuous  data  acquisition.  The  A/D  card  is 
externally  triggered  using  2-3  volt  TTL  type  signals  from  the  unsteady  device 
motor  drive.  When  the  unsteady  device  is  not  in  operation,  triggering  is 
provided  by  a  Hewlett-Packard  331  LA  Function  Generator  set  to  the  same 
frequency  as  that  used  with  the  unsteady  device  in  operation.  Sampling 
occurs  as  square  wave  voltage  changes  from  high  to  low. 

The  storage  capacity  of  the  A/D  buffer  is  64,000  bytes.  The  number  of 
imposed  unsteadiness  cycles  which  can  be  sampled  depends  on  the  desired 
number  of  data  points  per  cycle.  For  example,  if  100  points  per  cycle  are 


desired,  then  640  cycles  can  be  sampled.  In  this  study,  sampling  was  done  at 
100  points  per  cycle  for  all  cases,  where  cycle  refers  to  the  cycle  of  flow 
blockage  caused  by  the  vane  rotation  of  the  unsteady  device.  100  points  per 
cycle  were  chosen  so  that  there  would  be  sufficient  resolution  to  determine 
the  effect  of  unsteadiness  over  the  complete  phase  of  each  unsteadiness  cycle. 
e.  Datastorage 

The  Hewlett-Packard  microcomputer  is  used  to  acquire  the  data  from 
the  A/D  buffer,  using  a  data  acquisition  program  HOTWIREPAV,  described  in 
the  software  directory  (APPENDIX  A).  Raw  instantaneous  voltage  data  is 
acquired  and  stored  on  standard  3.5  inch  micro  diskettes  for  further 
processing.  Manual  inputs  to  the  program  include  bulk  flow  velocity,  flow 
blockage  frequency,  and  date  and  time. 

£  Data  Processing 

Processing  of  the  data  is  done  using  a  Hewlett-Packard 
microcomputer  using  the  programs  MAX_MIN  and  AVGPLOT.  This  software  is 
described  in  detail  in  the  software  directory  (APPENDIX  A).  Many  of  the 
curves  used  to  graphically  display  the  results  were  generated  using  a 
Macintosh  SE  personal  computer  and  the  graphics  program  CRICKETGRAPH. 
The  phase  averaging  technique  employed  by  the  program  MAX_MIN  is 
described  below. 

In  the  case  of  an  unsteady  flow  the  instantaneous  velocity  can  be  split 
into  three  components. 

U  *  0  +  U  +  u 

where  0  is  the  time  averaged  velocity,  u  is  the  imposed  phase  averaged 
velocity  and  Lf  is  the  instantaneous  fluctuating  component  u  is  the 


component  due  to  the  unsteadiness,  and  is  zero  for  steady  bulk  flow.  In  the 
present  study  u  and  D  are  combined  to  represent  the  phase  averaged  velocity 
(U\  Thus, 

U  *  u  +  if . 

u  is  determined  by  phase  averaging  U  over,  say  640  cycles  (ncycl),  at 
different  positions  of  the  phase  of  each  oscillation. 
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When  if  is  a  large  percentage  of  u  then  a  larger  number  of  cycles, 
sometimes  up  to  4640  cycles,  are  used  to  determine  u .  For  this,  multi-run 
data  acquisition  is  made,  with  the  vane  at  the  same  starting  position  for  each 
run.  Only  the  last  580  data  points  (due  to  sampling  procedure)  of  each  run  are 
used  in  the  cumulative  phase  averaging.  Once  u  is  determined,  if  is 
calculated  from  tf  *  U  -  a,  and  ((/)'  -  UJ  -  u)  .  Therefore,  the  phase 
averaged  root  mean  squared  velocity  is, 
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Time  averaged  velocity,  and  rms  velocity  \((/T  are  calculated  by 
averaging  over  all  the  points  ( 100  x  ncycl)  per  run  and  in  the  case  of  multi¬ 
runs,  100  x  ncycl  x  no.  of  multi-runs. 
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D.  FLOW  VISUALIZATION 

The  equipment  utilized  for  flow  visualization  studies  includes  a  smoke 
wire  for  generating  smoke  in  the  channel  and  video  equipment  to  photograph 
the  flow  in  the  channel.  These  are  described  in  detail  below. 

1.  Smoke  Wire  System 

The  smoke  wire  system  is  a  modified  Vertical  Smoke  Wire  Instrument 
made  by  Flow  Visualization  Systems  of  Bolingbrook,  IL.  The  system  in  its 
present  configuration  consists  of  a  0.127  mm  diameter  nichrome  wire  that  is 
stretched  across  the  width  of  the  channel  through  access  ports  in  the  channel 
side  walls.  The  wire  is  passed  through  these  access  ports  through  a  thin 
stainless  steel  tube  (hypodermic  needle)  which  allows  the  wire  to  move  freely. 
The  smoke  wire  is  connected  to  an  A/C-D/C  converter  and  this  in  turn  has 
voltage  supplied  by  a  Calrad  45-740  (0- 130  VAC)  variac  set  at  45  VAC.  The 
wire  is  energized  by  pulsing  the  "Momentary  On*  switch  of  the  A/C-D/C 
converter.  Since  the  wire  acts  as  a  resistor  and  heats  up  when  a  current  is 
passed  it  becomes  more  elastic  and  has  a  tendency  to  sag  in  the  channel.  To 
prevent  this  from  affecting  the  flow  visualization,  a  tensioning  device  has  been 
installed  on  one  end  of  the  wire  to  ensure  that  the  wire  is  held  taut  in  the 
channel  when  it  is  energized.  This  tension  device  consists  of  two  pulleys 


which  support  the  wire  and  a  third  pulley  with  a  weight  attached  between  the 
two  support  pulleys.  The  wire  is  coated  with  a  paraffin  based  oil  which  is 
applied  by  means  of  a  thin  rod  with  an  absorbent  piece  of  felt  at  the  tip.  When 
applied  in  this  manner,  the  wire  is  coated  evenly  with  small  droplets  of  oiL 
Energizing  the  wire  produces  a  thin  even  sheet  of  smoke  until  all  the  oil  on 
the  wire  bums.  The  smoke  from  the  wire  lasts  for  about  4-8  seconds 
depending  on  the  flow  rate  in  the  channel 

The  wire  is  placed  horizontally  normal  to  the  flow  76.8  cm  upstream 
from  the  end  of  the  test  section.  For  the  purposes  of  this  study,  the  wire  was 
located  at  a  y/d  value  of  0.84  (d  being  the  width  of  the  channel  and  y  the 
height  from  the  bottom  of  the  channel). 

2.  Video  Equipment 

The  video  equipment  used  for  flow  visualization  is  the  same  as  that 
described  by  Longest  [Ref.  9[.  Briefly,  a  high  speed  video  camera  is  mounted 
on  a  tripod  above  the  channel  and  focused  on  the  area  of  interest  The  focal 
plane  from  the  camera  lens  to  the  channel  is  surrounded  by  an  enclosure  to 
exclude  all  extraneous  light  and  thus  decreasing  the  reflections  which 
interfere  with  the  flow  visualization  The  bottom  of  the  channel  was  blacked 
out  and  lighting  was  provided  by  two  600  Watt  floodlights  placed  1  m  from  the 
sidewall  on  either  side  of  the  channel  at  the  same  height  as  the  channel 
Still  photographs  were  made  by  taking  pictures  of  the  aforementioned 
video  when  viewed  on  a  monitor.  A  Nikon  F-3  SLR  and  Kodak;Tri-Xpan  film 
were  used  to  take  the  pictures.  The  shutter  speed  was  one-eighth  of  a  second 
and  the  focal  aperture  was  set  at  f  5.6. 


HL  EXPERIMENTAL  PROCEDURES 


A.  ORIFICE  AP  CALIBRATION 

The  orifice  differential  pressure  calibration  for  bulk  velocity  flow 
measurements  was  performed  by  Longest  [Ref  9].  Results  of  the  calibration  are 
shown  in  Figure  4. 

B.  HOT-WIRE  CALIBRATION 

The  hot-wire  was  calibrated  in  the  freestream  of  the  wind  tunnel,  also 
located  in  the  laboratories  of  the  Mechanical  Engineering  Department  of  the 
Naval  Postgraduate  SchooL 

The  hot-wire  probe  is  mounted  normal  to  the  flow  in  the  center  of  the 
channel  so  that  it  is  in  the  freestream.  The  hot-wire  is  connected  to  the  same 
hot-wire  bridge  previously  discussed  Output  from  the  hot-wire  bridge  is  read 
on  a  digital  multimeter.  Freestream  velocity  in  the  wind  tunnel  is  measured 
utilizing  a  Kiel  pressure  probe,  a  wall  static  pressure  tap,  and  a  Validyne  digital 
manometer.  Voltage  and  differential  pressure  readings  are  taken  for  a  range 
of  pressure  drops  corresponding  to  velocities  of  1.0  to  4.0  m/s.  The  computer 
program  HWCAL  is  run  on  the  Naval  Postgraduate  School's  IBM  main  frame 
computer  to  calculate  all  calibration  constants.  Following  the  calibratioa  the 
hot-wire  is  installed  in  the  straight  channel. 


C.  CHANNEL  VALIDATION 
1.  Channel  Plow  Teats 

Assembly  and  installation  of  the  straight  channel  in  the  Mechanical 
Engineering  Laboratory  of  the  Naval  Postgraduate  School,  along  with  the  initial 
flow  visualization  and  measurement  results  are  detailed  by  Longest  [Ref.  9). 
Prior  to  commencing  flow  measurements  for  this  study  a  check  of  the  channel 
flow  quality  is  performed  as  follows.  Smoke  is  generated  using  the  smoke 
generator  described  by  Longest  [Ref.  9)  and  then  drfc*/n  into  the  channel  at  a 
speed  of  approximately  20  m/s.  The  first  few  times  this  was  performed, 
numerous  leaks  from  the  sidewalls  of  the  channel  were  observed.  These  leaks 
are  eliminated  using  removable  spring  clamps  on  both  sides  of  the  channel, 
along  its  entire  length.  The  clamps  are  left  in  position  while  the  channel  is  in 
operation. 

2  Laminar  Average  Velocity  Profile Verification 

Velocity  profile  verification  is  performed  at  a  Reynolds  number  of  1 103 
(bulk  velocity  of  1.32  m/s)  when  the  flow  is  laminar.  Instantaneous  velocity 
measurements  were  taken  at  eight  y/d  positions  ranging  from  0.2  to  0.9,  at  an 
x/d  of  300.0.  This  data  was  taken  using  the  HOTWIREPAV  program,  with  no 
imposed  unsteadiness.  MAX_M1N  was  used  to  process  the  data  and  the  time 
averaged  local  velocities  were  normalized  with  respect  to  the  maximum  value 
which  occurred  at  the  mid-channel  location  (y/d  =  0.5).  A  second  order 
polynomial  curve  was  fit  to  the  data  and  is  shown  in  Figure  5.  The  curve 
demonstrates  that  the  flow  is  fully  developed  at  this  location  and  follows  the 
expected  parabolic  behavior. 


0.  FLOW  MEASUREMENT 

L  Reynolds  Number  Survey 

Velocity  surveys  are  conducted  for  Reynolds  numbers  ranging  from  1 100 
to  2700  with  the  hot-wire  probe  at  a  y/d  position  of  0.85,  and  a  zjd.  of  -7.7. 
This  survey  is  done  to  determine  how  the  imposed  unsteadiness  affects  the 
transition  process  as  Reynolds  number  varies.  Data  at  each  Reynolds  number 
are  taken  by  recording  results  with  no  imposed  unsteadiness  followed  by 
results  with  imposed  unsteadiness.  Data  sets  of  64,000  points  are  taken  using 
only  the  last  58,000  for  data  processing,  as  previously  described.  Average 
velocity  and  rms  velocity  values  from  such  survey  are  used  to  generate  Figures 
6  and  7.  These  experiments  were  repeated  at  several  Reynolds  numbers  to 
ensure  consistent  and  repeatable  data  trends. 

2.  Multi-run  Averaging 

Multi-run  phase  averaging  using  up  to  4640  cycles  of  flow  unsteadiness 
were  conducted  at  several  Reynolds  numbers.  The  purpose  of  this  multi-run 
averaging  is  to  show  that  the  survey  data  is  statistically  stable  for  higher 
Reynolds  numbers  where  fluctuation  levels  are  significant  with  respect  to  the 
local  mean  velocity.  Results  are  presented  in  Figures  8-35. 

3.  Velocity  Profiles 

Velocity  profiles  with  and  without  imposed  unsteadiness  are  taken  at  the 
same  Reynolds  numbers  as  used  for  the  multi-run  averaging.  The  purpose  of 
these  profiles  is  to  record  flow  development  as  the  Reynolds  number  is 
changed  and  to  study  the  effect  of  imposed  unsteadiness  on  average  and  rms 
velocities  as  a  function  of  the  Reynolds  number  and  y/d.  Data  acquisition  and 


processing  is  performed  as  described  in  the  Reynolds  number  survey.  Results 
are  provided  in  Figures  36-147. 

E.  FLOWVISUALIZATION 

Smoke  patterns  were  recorded  on  film  using  the  video  camera  at  Reynolds 
numbers  ranging  from  1140  to  2240  with  1  Hz  imposed  unsteadiness.  Video 
sequences  are  taken  at  different  Reynolds  numbers  spaced  apart  by 
approximately  50  units.  The  purpose  of  this  flow  visualization  is  to  record 
transition  events  and  denote  any  variations  caused  by  the  imposed 
unsteadiness.  Still  photographs  are  shown  in  Figures  148-158. 


IV.  RESULTS  AND  DISCUSSION 


The  coordinate  system  used  in  this  study  is  a  right  handed  coordinate 
system.  The  origin  is  at  the  spanwise  centerline,  at  the  bottom  of  the  channel 
and  at  the  beginning  of  the  test  section.  The  x-axis  extends  positively  in  the 
streamwise  direction  to  the  end  of  the  test  sectioa  The  y-axis  extends 
vertically  in  the  positive  direction  from  the  bottom  of  the  channel  to  the  top. 
The  z-axis  extends  in  the  spanwise  directioa  Positive  values  of  z  are  on  the 
right  hand  side  of  the  channel  when  looking  downstream.  Flow 
measurements  were  taken  at  x/d  =  300,  z/d  =  -7.7  and  y/d  values  of  0.2  to  OS, 
where  d  is  the  full  channel  height  Flow  visualization  was  done  with  the 
smoke  wire  positioned  at  a  y/d  value  of  0.85. 

The  Reynolds  number.  Re,  is  based  on  the  bulk  flow  velocity.  Vbaik,  and  the 
channel  height  and  is  defined  by  the  equatioa 
Re  =  vboii:  •  d/v 

where  v  is  the  kinematic  viscosity.  The  Strouhal  number,  Str,  is  defined  by 
the  equatioa 

Str  =  2a  *  fee  *  d/vyxjk 

where  ft*  is  the  frequency  of  the  impc.ed  unsteadiness.  This  study  was 
conducted  at  Reynolds  numbers  from  1103  to  2715.  With  1  Hz  imposed 
unsteadiness,  Strouhal  numbers  range  from  0.0605  to  0.0246. 


A.  CHANNEL  FLOW  QUALIFICATION 

For  fully  developed  channel  flow  the  theoretical  velocity  relationship  is 
given  by  Schlichting  (Ref  ll:p.  280-28 1|  as: 
u/uomx  =  40(  y/d  -  (y/d)2] 

where  u  is  the  local  velocity  and  Um  is  the  centerline  (y/d  «  0.5)  velocity. 
Figure  5  shows  the  velocity  profile  taken  at  x/d  ■  300  at  a  Re  of  1 103.  The 
velocity  profile  (normalized  with  respect  to  Un*x)  is  parabolic  in  shape  and  is 
in  good  agreement  with  the  theoretical  relationship.  With  the  bulk  flow 
velocity  of  1.32  m/s  (measured  from  the  orifice  pressure  drop)  the  maximum 
centerline  velocity  of  1.774  m/s  is  observed,  giving  a  ratio  of  vw,iy/ut»»  of 
0.774. 

B.  EFFECT  OF  REYNOLDS  NUMBER  ON  LOCAL  MEAN  VELOCITY  AT 
y/d  =  0.85  DUE  TO  1  HZ  IMPOSED  UNSTEADINESS 

Figure  6  shows  the  variation  of  the  mean  velocity  with  Reynolds  number 
with  and  without  imposed  unsteadiness.  Flow  measurements  are  for  y/d  = 
0.85.  Bulk  velocity  is  also  shown  in  the  plot  The  mean  velocity  trace  shows 
some  non-linear  behavior  for  Reynolds  numbers  from  1900  to  2200  and  near  a 
Reynolds  number  of  2700.  This  behavior  is  believed  to  be  caused  by  the 
transitional  behavior  of  the  flow.  The  1  Hz  imposed  unsteadiness  has  no  effect 
on  the  mean  velocity. 


C.  EFFECT  OF  REYNOLDS  NUMBER  ON  RMS  VELOCITY  AT  y/d  =  0.85 

Figure  7  shows  normalized  longitudinal  velocity  fluctuations  with  and 
without  imposed  unsteadiness  for  Reynolds  numbers  from  1 100  to  2700. 

The  time  averaged  rms  velocity,  yfu*  ,  is  normalized  with  respect  to  the  time 
averaged  velocity,  V  ,  and  is  plotted  versus  the  Reynolds  number.  For 
Reynolds  numbers  of  1 100  to  1900,  the  imposed  unsteadiness  appears  to  have 
no  effect  on  the  rms  velocity.  For  Reynolds  numbers  from  1900  to  2025,  the 
1  Hz  imposed  unsteadiness  causes  a  significant  increase  in  the  rms  velocity. 
This  increase  amounts  to  a  maximum  at  a  Re  «  1980  and  is  7.06  %.  At 
Reynolds  numbers  of  2025  to  2300  the  effect  of  the  imposed  unsteadiness  is 
reversed.  In  this  range  the  imposed  unsteadiness  reduces  the  normalized  rms 
velocity.  The  maximum  amount  of  decrease  occurs  at  Re  *  2172  and 
corresponds  to  a  drop  of  7.3  %.  At  higher  Reynolds  numbers,  the  imposed 
unsteadiness  appears  to  have  no  effect  on  the  rms  velocity.  The  figure  also 
seems  to  demonstrate  that  transition  begins  near  Re  =  1200,  and  that  the  flow 
becomes  fully  turbulent  at  Re  =  2400. 

D.  EFFECT  OF  REYNOLDS  NUMBER  ON  PHASE  AVERAGED  VELOCITY 

Figure  8  through  Figure  35  show  the  phase  averaged  velocity  traces  at  five 
Reynolds  numbers  (1153,1826,  2005,  2085,  and  2640).  Figures  8-17  are  the 
phase  averaged  plots  with  no  imposed  unsteadiness.  Of  these,  figures  (  8,  10, 
12,  14,  and  16)  plot  phase  averaged  velocity  b(in  this  case  u  -  U  since  0  = 

0)  versus  phase  angle.  Figures  9,  11,  13,  15,  and  17  plot  the  phase  averaged 
rms  velocity,  y[tf*  t  normalized  with  respect  to  the  temporal  mean,  U ,  versus 
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phase  angle.  For  these  cases,  phase  averaging  is  done  over  580  cycles  of 
imposed  unsteadiness  with  100  data  points  per  cycle. 

Figures  18-35  present  results  of  flow  measurements  with  1  Hz  imposed 
unsteadiness.  These  figures  show  that  the  unsteady  device  produces 
deterministic  and  periodic  unsteadiness,  thereby  validating  its  behavior.  At  Re 
=  1 153  the  phase  averaging  is  given  for  580  cycles  of  imposed  unsteadiness. 
For  Re  =  1826,  2005  and  2085  eight  runs  of  580  cycles  each  (for  a  total  of 
4640  cycles)  are  averaged.  For  Re  *  2640  three  runs  of  580  cycles  (total  of 
1740  cycles)  are  averaged.  With  1  Hz  imposed  unsteadiness,  statistically 
stable  phase  averaging  is  achieved  using  580  cycles  when  the  flow  is  laminar 
(Re  =1153  for  the  plots  presented).  However,  for  the  transitional  Reynolds 
numbers.  Re  =  1826,  2005  and  2085,  eight  sets  of  580  cycles  are  necessary  to 
achieve  statistically  stable  phase-averaged  results.  For  the  turbulent  case.  Re  = 
2640,  statistically  stable  results  are  achieved  after  phase  averaging  three  sets 
of  580  cycles  of  imposed  unsteadiness.  The  effect  of  multi-run  averaging  is 
shown  in  Figures  20-35.  Figure  20  shows  u  after  averaging  using  one  set  of 
580  cycles.  Figure  22  shows  u  after  averaging  using  eight  sets  of  580  cycles, 
and  is  clearly  much  smoother  indicating  its  statistical  stability.  Figures  21  and 
23  show  the  same  effect  on  the  rms  velocity.  Figures  24-27  show  the  same 
results  of  multi-run  averaging  for  Re  =  2005,  Figures  28-31  for  Re  =  2085  and 
Figures  32-35  for  Re  =  2640. 

Normalized  rms  velocity  values  in  these  figures  are  fairly  constant  as  the 
phase  angle  changes  at  all  Reynolds  numbers,  with  the  exception  of  Re  = 

2085.  These  data  are  shown  in  Figure  3 1.  Here,  the  rms  velocity  increases 
for  the  phase  angles  of  250-350.  This  increase  occurs  at  the  same  phase 


angles  where  u  is  larger  than  stirrounding  values.  The  increase  in  rms 
velocity  values  may  occur  during  decelerating  portion  of  the  u  cycle,  followed 
by  a  phase  shift  relative  to  the  6/  trace.  These  rms  velocity  variations  occur 
over  the  same  Reynolds  number  range  that  imposed  unsteadiness  decreases 
the  rms  velocity  relative  to  the  flow  without  imposed  unsteadiness. 

EL  MEAN  VELOCITY  PROFILES  AND  PHASE  AVERAGED  VELOCITIES 

Figures  36-42  show  mean  velocity  profiles  at  x/d  =  300  and  z/d  =  -7.7. 
y/d  positions  vary  from  0.2  to  0.9.  Data  is  not  given  for  y/d  of  0.1  due  to 
physical  constraints  imposed  by  the  geometry  of  the  hot-wire  probe. 

Figure  36  shows  the  profile  development  with  Reynolds  number  with  no 
imposed  unsteadiness.  Profiles  are  shown  for  Re  =  1153,  1826,  2026,  2085 
and  2640.  With  laminar  flow,  profiles  are  parabolic.  Profiles  are  somewhat 
contorted  at  transitional  Reynolds  numbers  (Re  *  1826,  2026  and  2085).  At 
Re  =  2640,  the  profile  shows  a  flattened  behavior  expected  for  fully  turbulent 
flow. 

Figure  37  shows  the  profile  development  for  the  same  Reynolds  number 
with  1  Hz  imposed  unsteadiness.  The  imposed  unsteadiness  nominally  has  no 
effect  on  the  mean  velocity  profiles.  This  is  also  shown  in  Figures  38-42 
where  mean  velocity  normalized  with  respect  to  the  maximum  observed  mean 
velocity  with  and  without  unsteadiness  are  plotted  for  each  Reynolds  number, 
both  with  and  without  imposed  unsteadiness. 

Figure  43  shows  the  variation  of  the  peak  to  peak  magnitude  of  the  phase 
averaged  mean  velocity  as  a  function  of  y/d  On  the  plot,  the  peak  to  peak 
magnitude  has  been  normalized  with  respect  to  u .  From  the  figure,  the  peak 


to  peak  magnitude  increases  when  the  flow  is  near  the  walL  At  the  two 
transitional  Reynolds  numbers  (Re  -  2026  and  2085),  increases  in  magnitude 
near  the  wall  are  significantly  higher  than  for  other  cases.  These  Reynolds 
numbers  correspond  to  the  situation  in  which  the  imposed  unsteadiness 
decreases  the  rms  velocity  relative  to  flow  with  no  imposed  unsteadiness. 

Figures  44-91  show  phase  averaged  data  corresponding  to  data  just 
discussed  These  data  indicate  that  the  unsteady  device  is  successful  in 
producing  deterministic  periodic  unsteadiness  at  all  positions  where  flow  is 
measured,  i.e.  from  y/d  =  0.2  to  0.9. 

F.  RMS  VELOCITY  PROFILES 

Figures  92-98  show  the  development  of  the  rms  velocity  profiles  as  the 
Reynolds  number.  Again,  the  symmetry  of  the  profiles  is  distorted  due  to  the 
fact  that  no  measurements  are  available  at  y/d  *  0.1. 

Figure  92  presents  rms  velocity  profiles  for  the  five  Reynolds  numbers 
with  no  imposed  unsteadiness.  Figure  93  shows  the  profile  development  with 
imposed  unsteadiness  at  1  Hz.  Figures  94-98  are  normalized  rms  velocity 
profile  plots  with  and  without  imposed  unsteadiness  for  each  Reynolds 
number. 

Figures  92  and  93  show  that  the  imposed  unsteadiness  has  no  effect  on 
the  rms  velocity  profiles  for  the  laminar  and  turbulent  Reynolds  numbers  (Re 
«  1153,  1826  and  2640).  For  the  turbulent  case,  the  rms  velocities  are  much 
higher  near  the  wall  For  the  two  transitional  Reynolds  numbers,  some 
different  behavior  is  seen.  With  no  imposed  unsteadiness,  the  profiles  show 
higher  rms  velocity  magnitudes  near  the  middle  of  the  channel.  With  imposed 
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unsteadiness,  the  profiles  at  Re  *  2026  and  Re  =  2085  differ  sharply  and 
resemble  the  turbulent  profile. 

Figures  99-146  present  phase  averaged  data  used  to  produce  the  plots 
described  above.  Rms  values  are  constant  as  a  function  of  phase  angle  both 
with  and  witnout  imposed  unsteadiness  except  for  the  case  with  imposed 
unsteadiness  at  Re  *=  2085.  In  Figures  116,  118,  120,  12Z  124,  126,  128, 
and  130  rms  velocities  vary  with  phase  angle.  This  was  mentioned  earlier  for 
y/d  »  0.85,  however,  here  it  is  evident  that  this  behavior  occurs  for  0.2  <  y/d  < 
0.9 .  Figure  147  shows  the  variation  of  rms  peak  to  peak  magnitude  as  a 
function  of  y/d  for  Re  a  2085.  At  this  Reynolds  number,  the  peak  to  peak 
magnitude  is  largest  near  the  wall  with  a  secondary  peak  near  the  channel 
centerline. 

G.  FLOW  VISUALIZATION 

Figures  148-158  show  flow  visualization  results  obtained  with  the  smoke 
wire  at  y/d  of  0.85.  The  camera  was  set  up  to  take  videos  on  the  left  side  of 
the  channel  in  the  same  area  that  the  hot-wire  was  located  (x/d  =  300).  In 
the  figures,  flow  is  from  the  bottom  of  the  picture  to  the  top.  The  left  side  of 
the  picture  is  near  the  left  channel  wall,  and  the  right  hand  side  is  located 
near  the  center  of  the  channel.  The  dark  horizontal  stripe  emanating  from 
the  left  sides  of  the  pictures  indicates  the  location  of  the  hot-wire  probe. 

There  are  also  several  dark  shadows  which  extend  diagonally  upward  from  the 
left  side  of  the  channel.  These  are  due  to  shadows  cast  by  part  of  the  support 
structure  of  the  channel  and  should  not  be  confused  with  effects  of  the  flow  on 
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the  smoke  patterns  .  The  smoke  generated  by  the  wire  appears  as  white  lines 
or  ribbons  and  the  dark  areas  indicate  absence  of  smoke. 


Figure  148  shows  flow  at  Re  =  1 136  with  1  Hz  imposed  unsteadiness.  In 
this  picture,  the  flow  is  fully  laminar  as  shown  by  the  straight  regular  parallel 
smoke  lines.  Figure  149  is  at  Re  -  1 153  also  with  imposed  unsteadiness. 

Here  flow  is  again  laminar,  but  there  are  also  TS  (Tollmien-Schlichting)  waves 
in  addition  to  some  small  streamwise  spirals  due  to  secondary  vortical  flows. 
These  are  hereafter  referred  to  as  ribbons.  Figure  150  at  Re  *  1278  again 
shows  laminar  flow  with  TS  waves.  Figure  15 1  at  Re  *  1328  also  shows 
laminar  flow  with  TS  waves.  Figures  152  and  153  are  at  Re  -  1395,  and  show 
TS  waves  and  ribbons  which  are  larger  than  previously  seen  and  which  begin 
to  spread  out  diagonally.  Figures  154  and  155  are  at  Re  =  1428  again  show 
the  same  ribbon  structure,  but  here  they  have  become  larger  and  seem  to 
develop  into  turbulent  spots.  These  spots  are  defined  by  the  large  circular 
areas  in  the  center  of  the  ribbon  break-up  regions.  Figures  156  and  157  are 
for  Re  =  1586.  These  pictures  show  the  decay  of  ribbon  structures  as  well  as 
larger  and  more  frequent  turbulent  spots.  Figure  158  shows  a  time  sequence 
of  six  pictures  taken  approximately  one  second  apart  for  Re  =  1428.  In  this 
sequence  we  see  laminar  flow  with  TS  waves  is  present  in  addition  to  ribbons, 
ribbons  developing  into  turbulent  spots,  followed  by  more  ribbons  and  laminar 
flow  with  TS  waves.  The  principal  difference  between  the  flows  with  and 
without  imposed  unsteadiness  is  the  reappearance  of  such  TS  waves  at  higher 
Reynolds  numbers  with  imposed  unsteadiness.  In  this  situatioa  the  TS  waves 
generally  reappear  just  downstream  of  turbulent  spots. 
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V.  CONCLUSIONS 


The  1  Hz  imposed  bulk  flow  unsteadiness  has  no  effect  on  the  local  mean 
velocity  at  all  y/d  values  for  the  range  of  Reynolds  numbers  1 100  to  2700. 

The  imposed  unsteadiness  causes  both  destabilization  and  stabilization  of 
the  time  averaged  longitudinal  velocity  fluctuations.  The  flow  is  destabilized  by 
the  imposed  unsteadiness  for  Reynolds  numbers  1900  to  2025  and  Strouhal 
numbers  from  0.0351  to  0.0329.  Stabilization  of  the  flow  by  the  unsteady 
device  occurs  for  the  Reynolds  numbers  2025  to  2300  and  for  Strouhal 
numbers  from  0.0329  to  0.0290. 

Imposed  unsteadiness  affects  the  shape  of  the  time  averaged  longitudinal 
velocity  profile  at  the  transitional  Reynolds  numbers.  Re  *  2026  and  Re  * 

2085. 

The  imposed  unsteadiness  also  causes  the  time  averaged  longitudinal 
velocity  fluctuations  to  vary  as  a  function  of  phase  angle  at  a  Reynolds  number 
of  2085  at  y/d  positions  of  0.2  to  0.9.  The  maxima  of  the  time  averaged 
longitudinal  velocity  fluctuations  and  the  maxima  of  the  phase  averaged 
velocity  occur  for  the  same  phase  angles. 

Transitional  flow  structures  which  have  been  predicted  by  numerical 
simulations  and  observed  in  other  experimental  work  are  also  observed  in  this 
study.  As  Reynolds  number  increases  laminar  flow,  Tollmien-Schlichting 
waves,  vortical  ribbon-like  structures  and  then  turbulent  spots  are  observed. 
These  structures  are  observed  to  be  present  at  many  different  Reynolds 
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numbers.  Time  sequence  photography  taken  at  particular  Reynolds  numbers 
during  transition  shows  TS  waves,  vortical  ribbons,  turbulent  spots,  more 
vortical  ribbons  and  TS  waves.  The  growth  and  break-up  of  the  ribbon-like 
vortical  structures  is  followed  by  turbulent  spots.  As  the  Reynolds  number  is 
increased  the  vortical  structures  grow  in  both  the  spanwise  and  streamwise 
direction.  The  frequency  and  size  of  the  turbulent  spots  also  increases  as 
Reynolds  number  increases  until  the  flow  becomes  fully  turbulent  The 
principal  difference  between  the  flows  with  and  without  imposed 
unsteadiness  is  the  reappearance  of  the  TS  waves  at  higher  Reynolds  numbers 
with  imposed  unsteadiness.  In  this  situation,  the  TS  waves  generally  reappear 
just  downstream  of  turbulent  spots. 


APPENDIX  A 


SOFTWARE  DIRECTORY 

1.  HWCAL:  This  program  determines  the  constants  for  the  King's  Law 
calibration  of  the  hot-wire.  The  program  also  provides  a  polynomial  fit  of  the 
calibration  data 

2.  HOTWlKEPAV  :This  program  is  used  to  read  the  data  stored  in  the  A/D 
buffer  of  the  high  speed  data  acquisition  system  and  stores  the  information  on 
micro  diskettes.  Manual  inputs  are:  triggering  frequency,  hot-wire  DC  voltage 
(ungained),  oscillation  frequency  (flow  blockage),  bulk  velocity  and  date  and 
time  of  run. 

3.  MAX-MIN :  This  program  calculates  instantaneous  and  phase  averaged 
velocities.  Initially  a  look  up  table  is  created.  Here,  effective  velocities  are 
calculated  from  the  effective  voltage  values  and  stored  for  follow  on 
calculations.  The  hot-wire  calibration  constants  obtained  from  HWCAL,  and 
the  amplifier  gain  are  incorporated  into  these  calculations.  The  velocity 
calibration  is  given  by  the  equation: 

Uefr«  k(E*r2-Eo2)i/N 

where  k  is  the  proportionality  constant,  is  the  effective  voltage,  and  E0  is 
the  reference  voltage  at  no  flow.  N  is  a  constant  value  of  0.45  for  moderate 
Reynolds  numbers.  Once  the  look  up  table  is  created  the  program  reads  the 
instantaneous  voltage  values  from  the  data  file  and  converts  them  to 


instantaneous  velocities.  At  this  point  a  plot  of  the  instantaneous  velocity 
versus  time  can  be  generated  for  any  of  580  cycles  (The  first  60  cycles  are 
discarded  to  allow  for  flow  stabilization  as  the  unsteady  device  is  started) 

Next  the  program  phase  averages  the  580  cycles,  and  velocity  versus  phase 
angle  (of  the  flow  blockage)  plots  are  available.  Two  plots  are  available  from 
the  averaged  values;  u  versus  phase  angle  and  divided  by  U,  versus 
phase  angle,  where  ti  is  the  phase  averaged  velocity,  VP  is  the  phase 
averaged  root  mean  squared  velocity,  and  U  is  the  average  velocity.  In  the  case 
where  there  is  no  imposed  unsteadiness,  the  phase  averaged  velocity,  u,  is 
equal  to  the  time  averaged  velocity ,  u  since  there  is  no  superimposed  phase 
average  velocity  component,  u. 

The  program  can  also  be  used  to  process  a  set  of  data  when  multi-run 
averaging  is  desired  In  this  case  the  averaged  values  are  written  to  and  stored 
in  the  computer  hard  disk  memory.  Once  each  of  the  desired  data  sets  has 
been  averaged  and  stored  in  the  computer  memory. 

4.  AVGPLOT :  The  program  AVGPLQT  is  used  to  average  any  number  of  data 
sets.  A  maximum  of  eight  sets  of  data  (or  4,640  cycles)  were  averaged  for  this 
study.  The  program  AVGPLOT  provides  the  same  graphs  as  those  available 
from  MAX_MIN.  Phase  averaged  data  sets  are  copied  to  micro  diskettes  prior 
to  purging  them  from  the  computers  memory. 
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FIGURES 


SCHEMATIC  OF  FLOW  MEASUREMENT  EQUIPMENT 


Figure  3. 


36 


Figure  4, 
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FullyDeveloped  Flow  Validation 
at  Re  No.  =  1153 
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Figure  8. 


41 


PHRSE  ANGLE  (dog. 


N2J  130999.  1938 


J  a 

Lu  f~ 

>  cn 

U) 

£  N 

K I 


-  a 

in  u. 


Figure  9. 


rrrrrirrrrjrTrrrrTTT;  in  rpm jTmj rrrrr 


L.uj-Lujj-1-U uLlj  ixLuxi! i.;.M.Ia_u.i.Uaiiiu  ixlxuiL 

I  *  80  *  90  *  MT  20  *  0 

~!3A  -DAy.'iyH.n 


PHASE  ANGLE  (deg. 


1  ? 09  8  9 . 2  04  ? 


W 


\ 

CD 

£ 

OJ 

m 

CD 

in 

-r-H 

CD  Q 

^ — i 

• 

# 

# 

2 

OJ 

2 

ZJ 

OH 


_J  •  >- 
LJ  £  Ld 
>  h  CK 
0) 

U) 

2 

£  N  1A 

I  \ 
Q  £ 
«  in 

\  *  CD 

£  0  — 
N  Li  . 
Q  Of  oj 
CD  h_ 

•  • 

^  •  -J 
V  LJ 

j  m  > 

LJ  O 
>  Y 
J 
D 
PQ 


0 

> 

cr 


iJ4  I  1 


(S/ui)  >IyHn 


Figure  10. 


RLJN2L  150989.2201 


S3 

CO 

CO 


Figure  12. 


45 


PHASE  ANGLE  (deg. 


0 

\ 

£  in 

«-<  Q 

ID  Q 

SQOJ 
cn 


#  # 


Q  J  •  ■ 

oj  y  (v  > 

ru  ! —  Ld 
(n  d 

m  in 
00  Z 
cd  n:  n 
(3  I  0 
in  Q  \ 
0  £ 
\  •  ^ 
£  0  • 
_l  ru  Ld  oj 

OJ  CD  q: 

Z  QIl  * 
d  ■  j 
&  cm  •  hi 

u  > 


j  tn 

u  o  y 

> 

J 

D 

• 

ra 

(J 

> 

d 

a 


Figure  13. 


46 


PHRSE  RNGLE  (deg. 


RUN3L  150989.2354 


\  in  L 

E  CO 

GO  0 

in  oj 

in  0  I 

cn 

•  # 

#  H 


J  •  > 
LxJ  £K  LxJ 

>  h-  q: 

V) 

U) 

T 

CK  N  CO 
I  \ 
0  E 


J- 

h 

V 


co  in 
\  •  0) 
EO'f 
OJ  Ld  • 

0  n  oj 

OJ  L_ 


OJ  •  J 
O  LJ 

j  tn  > 
u  o 

>  y 

j 

•  3 

CD  CQ 

> 

£E 


K 


rrrr,7n  ,'TP'rTrrnTT 


6  d 


■  ■  *  ■  t1  fe J 


Figure  16. 


iBHn 


49 


w 

\ 

\ 

£ 

w— < 

Q  ® 
'?• 


G 

CD 

ai 


2  01 
LI  O  v 
-> 


> 


(J 

> 

d 


* 

J 

D 

rc 


s 

co 


rrTTTprnjj  rypr  rrrrrrT  rp  m  pTTTp-rrrpTr  ?  jn7T|CO 

t  C  -1 

f- 

* — » 

(- 

t 


CD 

# 

# 

• 

r 

h 

i~ 

f- 

IT 

}  • 

> 

LI 

h 

00 

lj  a 

- 

'>  i_ 

r- 

■ 

✓  i 

01 

a 

f 

H- 

01 

01 

L 

00 

z 

►- 

Cl 

&  N 

w 

\ 

L 

(3 

2 

h 

<3 

Q 

\ 

E 

r 

r 

OJ 

w 

\  • 

CD 

u 

f- 

01 

£  O 

rH 

• 

1 

r 

"t  Li 

01 

▼— < 

i. 

2 

rv  L. 

i 

5 

• 

J 

a 

OJ  • 

i  1 

LJ 

L 

i- 

r 


Lu 

2' 


v 

CVl 

m 


-cd 

-  CO 
_  ai 


lai 

-iin 

-,ru 


JjCD 

4ai 

i 

1 

4 

12 

•4 

_f 

J? 


I00 

—4  gj 

u 

U 

Jcu 


j 

ju> 


lLu  u.i  j  GjIujuUuiIi  mlxi  u1j.jj  xinaai  i_LtxJ 

91*  2T  '  80  *  *0* 


tP 


13A  'SAU/ibH,"  SWi! 


Figure  17. 


50 


PHRSE  RNGLE  (deg. 


RUN2 I  13  0989.  1  9  2S 


jTiTqTrtrpTrrjTrn-jTrrq 


rTjTTTTj  n  rrpvrrjto 

4 

j 


\ 

\ 

E 

CD 

Q 

Q 


n 

oj  in 

03  - 
(\.  — 
in 
0 
•  # 


J  #  > 
U  hi 

>  •  a 

a 


tn  h* 

K 

I 

r  u) 

t 

a  ia 

\ 

\ 

L 

M  £ 

p 

W  I  CD 

*- 

i 

\h(I) 

t— 

E 

r 

rH  i  t~4 

r- 

m  O 

r~ 

u 

OJ  Ld  ■ 

„ 

03  a  j 

y- 

•  a  a 

- 

% 

L. 

** 

j  • 

i- 

hi  o  y. 

y 

>  G°>  J 

*— 

o  a 

•  m 

L 

(j 

L 

% 

L 

c 

i 

f- 

Lu  lLuj  d  uj  ili-iiU-ii  i  i  li-t  i  rfu-uJ-^-ulj-uila-iuJro 
8SE6  *  HEG*  ^926  ’  8126*  UI6*  S2I6' 


(S/wj  iyHn 


Figure  18. 


51 


PHASE  RNGLE  (dog 


RUN2 I  1 3  0989 . 1 92 S 

flVG  .  VEL  .9231  m-'s  RMS  VEL_  .00440m/s 

OSC.  FREQ.  1  Hz  STR .  #  .05782 


6} 


n 

in 


* 


y 

a 


M 

\ 

E 

CD 

ro 


Lu 


_ ! 

□ 

m 


Figure  19. 


iTryrrrry  ti  ttjttj'  rpri  i  jTTTTyrTrryn  rt  jtttt  yrnrntD 


r 

u 

r 

r 

r 

H 

L 

! 

h* 

t 

i 

r 

'r 


f- 

L. 

t" 

w 

L 


j- 


r 

h 


L 

L 


IV 
' — T  CVI 

4 


CD 

03 

ru 


;ru 

-Tin 

-,cu 


1 

4oj 


ia 

~~1CD 


I 

“M- 
—  ^r 
-4*~* 

J 

4 

lao 

~ia 


iW 


CD 

'In 


Llu  ilj-u-u  Lu  i-bxu  Liu  iiu-uiia  uluxilixuli-uaira 

I*  80*  90*  *0*  20-  0 

13A  *DAU.^iyH<n  swy 


52 


PHASE  RNGLE  (deg. 


PHRSE  RNGLE  (deg. 


(  2  M ,  1  RUNS)  PH.  RVE.  ,  170989 


(/> 

\ 

E  OJ 
in  in 
CD  co 
in  cn 
—  Q 


r  n  ryrTTrp-rm 


z  u 
□  > 
[ Y 


H  U 
J  > 
D  d 
Z 


“13A  *  DAb/ldH ,  n  SWa 


Figure  21. 


PHRSE  RNGLE  (deq. 


Figure  22. 


PHRSE  RNGLE  (deg. 


RUNS)  PH.  RVE. , 170989 


(A 

\ 

E  OJ 

n  in 
cd  CD 
-  CD 

•  0 


CD 

OJ 

CD 


# 


J 

LJ  # 

> 


> 

u 

a 


CO  a 
X  H 

a  in 


tfl  N 

\  I 
E  - 

CD  (D 

0  • 

-  CD  O 

X  •  Id 

OJ  -  CK 

~  L 

l 

Z  LJ  • 

□  >  U 

a.  in 

H  •  O 
H  (J 
J  > 

□  CE 

X 


w 

\ 

E 

in 

CD 


OJ 


J 

LJ 

> 

Y 

J 

□ 

m 


<s 


cn 

© 

T3 
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